Spin relaxation and band excitation of a dipolar BEC in 2D optical lattices 
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We observe interband transitions mediated by dipole-dipole interactions for an array of ID quan- 
tum gases of chromium atoms, trapped in a 2D optical lattice. Interband transitions occur when 
dipolar relaxation releases an energy larger than the lattice band gap. For symmetric lattice sites, 
and a magnetic field parallel to the lattice axis, we compare the measured dipolar relaxation rate 
with a Fermi Golden Rule calculation. Below a magnetic field threshold, we obtain an almost 
complete suppression of dipolar relaxation, leading to metastable ID gases in the highest Zeeman 
state. 

PACS numbers: 03.75.-b, 67.85.-d, 37.10.Jk 



Collisions and interactions between cold atoms trapped 
in optical lattices are subject to intense investigations. 
These studies are for example relevant for novel strongly 
correlated quantum phases [l[, and for metrology 
For most of the experiments to date, focus was made on 
contact interactions between atoms in the lattice vibra- 
tional ground state. Dipole-dipole interactions (DDIs), 
which have attracted considerable attention, especially 
since the recent production of Bose Eintein condensates 
(BECs) with highly magnetic chromium atoms |3|], in- 
troduce important new features in this context, because 
they can couple atoms in different lattice sites due to 
their long range character Here we show that DDIs 
can be used as well to couple different bands in the lattice, 
when magnetization changing collisions |5| release an en- 
ergy overcoming the gap between bands. In addition to 
their interest for quantum computing, higher lattice or- 
bitals may be used to study novel quantum phases, as 
suggested in [6]. 

Magnetization changing collisions in a cylindrically 
symmetric system lead to mechanical rotation in analogy 
to the seminal work of Einstein and de Haas (EdH) [7j , 



as suggested by [81-41 11] - Here, we study dipolar relaxation 
of spin S — 3 Cr BECs loaded in a strongly confining 2D 
optical lattice defining an ensemble of ID quantum gases 
in cylindrically symmetric sites. We observe a threshold 
as a function of the magnetic field, corresponding to the 
resonant exchange of a quantum of spin excitation and a 
quantum of rotational excitation in the lattice. Below the 
threshold, dipolar relaxation does not release enough ki- 
netic energy to induce the change in angular momentum, 
and we observe an almost complete suppression of dipo- 
lar relaxation for atoms in the highest energy Zeeman 
state (ms=+3); ID quantum gases in this state remain 
well below the degeneracy temperature 12| for tens of 
ms. Above threshold, dipolar relaxation occurs, with a 
transfer of population into excited bands of the lattice. 
While rotating states in each lattice site are then pre- 
dicted [ill ], and the observation of such rotation would 



be a proof of the EdH effect, these rotating states are 



not reached at our lattice depths due to fast tunneling 
between sites in excited band states. 

In our experiments, almost pure Chromium BECs in 
the absolute ground state ms = — 3 |l3|, with 20 000 
to 30 000 atoms, are loaded into a 2D optical lattice. 
The lattice is realized by two independent and orthogo- 
nal retro-reflected ID optical lattices, detuned from each 
other by 160 MHz. All beams derive from a 532 nm solid 
state single line laser. As the standing wave is turned 
on sufficiently slowly (within 15 ms), the BEC is loaded 
adiabatically in the lowest energy band of the 2D lattice, 
whose maximal depth in each direction is about 25 Er, 
where En = 2 m\i ^ s the recoil energy, with A = 532 nm, 
and m the mass of 52 Cr atoms. The energy difference in a 
lattice site between the fundamental and the first excited 
state, cji/27r=120 kHz, is much larger than the chemical 
potential, and the motion is frozen along two dimensions. 
The length of the harmonic oscillator ground state wave- 
function in the lattice sites = y^fi/muiL is about 40 
nm, and the peak density is raised up to 3 x 10 21 m~ 3 . 
As the tunneling time from one lattice site to the next is 
long (40 ms), we consider that the BEC is split into an 
array of a few hundred independent ID quantum gases. 

We then reduce the magnetic field B, so that the Lar- 
mor frequency ujq = g.jHsB/h (g,j is the Lande factor, 
He is the Bohr magneton) is adjusted close to ojl- The 
direction of B, u z , is set parallel to the direction of the ID 
tubes, ul («.e. perpendicular to the two standing waves 
axes). With a radio- frequency sweep, we transfer the 
atoms in the ms — 3 state: dipolar relaxation then be- 
comes energetically allowed. 

To characterize dipolar relaxation, we implement the 
following measurement procedure. We measure the num- 
ber of atoms in the different bands of the lattice by reduc- 
ing the lattice depth in 200 fis, which is slow compared 
to the timescale for band excitation in the lattice (10 (is), 
but fast compared to the timescale for thermalization (on 
the order of 1 ms) : this band mapping procedure [14j adi- 
abatically transfers the quasi-momentum distribution in 
the lattice into the real momentum distribution, which 
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FIG. 1: (Color online) Band mapping results above threshold. The 
false colored picture is an average of 4 absorption images of the BEC 
released from the lattice, (a) Integrated population profile along z: the 
first excited band (2 nd Brillouin Zone, BZ) in the y direction is popu- 
lated (in red), contrary to below threshold (in grey); noisy data points 
on the right (due to fringes in the absorption image) were suppressed, 
(b) Integrated profile along y showing a non Gaussian velocity distri- 
bution; lines arc results of a double Gaussian fit, yielding an effective 
temperature, (c) Sketch of the 2D lattice arrangement: the ID traps 
are along z, the imaging beam is sent along x. 
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FIG. 2: (Color online) Evidence for a threshold in dipolar relaxation: 
(a) v — 1 population as a function of Larmor frequency after 25 ms of 
dipolar relaxation for 25 E T \ (b) Effective temperature along the tubes 
as a function of Larmor frequency after 75 ms, for two different lattice 
depths (black circles: 12 E r : red squares: 25 E r ). The thick grey 
shaded area and the narrow red shaded area represent the respective 
widths of the second energy band. Lines arc guides for the eye. 



we image after a time of flight equal to 5 ms. The laser 
beam for absorption imaging being perpendicular to the 
tubes, we image the population of the different bands of 
one of the ID lattice (along y), and the velocity distribu- 
tion along the tubes (axis z), as shown in Fig. [TJ 

When B < B t h — fiVL/ '{9j^b), we observe almost 
no heating along the lattice tubes, and no population in 
higher bands. Typical results for magnetic fields above 
B t h are represented in Fig. [TJ We observe population in 
the first excited band (v — 1), and, within our signal to 
noise limit, no population in the second one (v = 2). In 
addition, we also observe a strong heating of the cloud 
along the z axis. Population in v — 1 and the effective 
(see below) temperature along the tubes are both repre- 
sented in Fig. [2] as a function of B. 

To explain these results, we recall that DDI between 



two atoms of magnetic moments /ii = gj^ssi and ji2 
9J/J-BS2 (si=(i,2) are spin operators), reads: 



V dd (r) 



fJ-o (ffj^a) si • s 2 - 3 (si ■ u r ) (s 2 ■ u r ) 
47r r 3 



(1) 



where f — ru r is the relative position of the two atoms, 
and /ito is the vacuum permeability. Vdd is invariant under 
simultaneous rotation of spin and space coordinates along 
any axis, so that the quantum number mj = ms + ttil 
associated to the projection of the total angular momen- 
tum S z +L z is conserved in dipolar collisions. In presence 
of B, when u z \\ ul, and for the specific case of cylindri- 
cally symmetric lattice sites, mj remains a good quan- 
tum number: dipolar relaxation results in the creation of 
rotating states of quantized orbital momentum. 

There are two distinct channels of dipolar relaxation 
when atoms are in ms 
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The first channel corre- 
sponds to Airig = — 1, the second channel to Amg = —2, 
and the corresponding release of kinetic energy is re- 
spectively Tiloq and 2haj a . Due to conservation of an- 
gular momentum, the first (second) channel associated 
to Ami = 1 (2) creates a rotating state of pairs (x + iy) 
((x + iy) 2 ), with x = x\ — x 2 and y = yi — yi the rel- 
ative coordinates. The corresponding energy costs are 
respectively Ulol and 2Tilol- As a consequence, both 
dipolar relaxation channels are expected to become en- 
ergetically possible for ujq > u>l, which explains the 
threshold observed in Fig. [2j Dipolar relaxation above 
threshold therefore involves two routes, which both pop- 
ulate excited vibrational states in the lattice. Chan- 
nel 1 populates the first excited lattice band, and, as 
(x + iy) 2 = (xi+iyi) 2 -2(xi+iyi)(x 2 + iy2) + {x2+iy2) 2 , 
channel 2 populates both the first and the second excited 
bands, and creates (motional) entangled states. 

As shown in Fig. [Hb), the threshold for dipolar re- 
laxation does depend on the lattice depth, and it corre- 
sponds to lo > lo l = 2ir x 120 (80) kHz for 25 (12) E R . 
Furthermore, the width of the observed threshold is com- 
parable to the width of the second excited band in the 
lattice (equal to 16 (43) kHz for 25 (12) E R ). As channel 
2 is the only one leading to population of the second ex- 
cited band we deduce that most of our signal comes from 
it. Finally, Fig. [2jb) also shows that raising the lattice 
depth leads to stronger heating above threshold. Such an 
increase cannot be solely interpreted as the effect of an 
increased density (the density only increases by only 40 
percent while dipolar relaxation roughly doubles from 12 
E r to 25 E r ), and will be related below to the reduction 
of the width of the second excited band. 

We interpret the strong heating along the tubes above 
threshold as the consequence of collisional de-excitation 
(vibrational quenching due to collisions with atoms in 
v — 0, see 16]). Due to conservation of orbital momen- 
tum, vibrational quenching is not possible for rotating 
states; however, simple calculations show that a singly 
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or doubly quantized vortex does not survive to a tunnel- 
ing event. As the tunneling time in the excited bands of 
the lattice is fast (1 ms in the first band, and 50 fis in the 
second band for 25 E r ), atoms rapidly lose their orbital 
momentum. A symmetry approach shows that v = 2 
atoms can undergo collisional de-excitation with atoms 
in v = (contrary to atoms in v = 1) [l6j ]: an atom in 
v = 2 can collide with an atom in v — 0, producing two 
atoms in v = 0; therefore population accumulates only in 
the first excited band, as observed in Fig. [2j As the pop- 
ulation in v = 1 rises, vibrational de-excitation between 
two colliding v = 1 atoms becomes efficient. Vibrational 
de-excitation from both bands therefore contributes to 
the observed heating; this heating is an indirect proof 
that atoms cease to rotate due to tunneling. 

An important feature of our experimental data is that 
the velocity distribution along the tubes rapidly becomes 
non-gaussian (see Fig. [TJ). Collisional de-excitation pro- 
duces pairs of back-to-back moving atoms at a relative 
momentum set by the lattice depth, so that the created 
velocity distribution is highly out of thermal equilibrium. 
It would be expected that it remains so in a pure ID ge- 
ometry with only contact interaction, because of the in- 
tegrability of such systems [ItJ • But in our situation the 
system is not purely ID, as the de-excited atoms have 
enough energy to populate higher bands in the lattice; 
and DDI may also come into play [3]. Despite this lack 
of integrability, the observed velocity distribution clearly 
departs from thermal equilibrium even after tens of colli- 
sions per particle. This intriguing lack of thermalization 
also has an impact on how we measure heating rates in 
our system: instead of fitting the velocity distribution 
along z by a simple gaussian, we evaluate its second mo- 
ment, hence defining an effective temperature. 

With the measured heating rate, and the measured 
population in v = 1, we deduce the total rate of ki- 
netic energy deposited in the system by dipolar relax- 
ation above threshold. We plot in Fig. [3] the population 
in v = 1 and the effective temperature as a function of 
time. In this experiment, the atoms are first promoted 
to m = 3 while B is below the threshold, and then B 
is raised slightly above the threshold (the threshold is 
reached at t = 25 ms in Fig. [3]). Below threshold, we 
observe no heating, and above threshold, we measure an 
energy increase rate of 215(±30) nK ms -1 . 

To account for these results and estimate a rate pa- 
rameter, we develop a model describing two particles 
in a Gaussian trap, assuming that dipolar relaxation 
can be described by a local density dependent rate (as 
in the 3D and 2D cases, see [IH). DDI at threshold 

couples a two-body state ^ cx exp (^—T^j exp t^?) 
describing a pair of atoms in the ground state of a 

lattice site (z = {-^f" » Po = (tBi)^) to 
rotating states which are excited not only perpendic- 
ularly but also along the lattice tubes. In the case 
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FIG. 3: (Color online) Dynamics of dipolar relaxation when the 
magnetic field is increased above threshold. Threshold is reached at 
t— 25 ms, (a) shows the effective temperature (see text), and (b) the 
measured proportion of atoms in the first excited band. Lines arc guides 
for the eye. 



of channel 2, the matrix element of DDI writes as 
^£exp(2i(9) [H, where d 2 = ^j^b) 2 . the ground 
state is only coupled to even vibrational states along 
the tubes, corresponding to wave functions ^2 



P 



exp(2i6>) exp (-^rj exp (-^) H 2n (^) (where H n 

stands for the n th Hermite polynomial). Similarly, 
channel 1 couples to odd vibrational states "J" cx 

pexp(*0) exp {-£r) exp H 2n +i (-|) . As an ex- 

ample, we calculate (for zq 3> po): 



V? = (*„ \V M \ *?) = 

V 2 ° = (H! \V M \ *§) = 



1/2 



S 3 / 2 d 2 



Po 



2 

Pq z Q z 



(27T) 



1/2 



Sd x 



P0Z0 



(2) 
(3) 



We see that V® -C V 2 °, as zo ^ Po : a t threshold, dipolar 
relaxation is therefore dominated by channel 2. This is 
consistent with the observed fact that the width of the 
threshold is comparable to the width of the second ex- 
cited band of the lattice. 

We stress that for our values of the lattice depth, V§ 
(resp. Vi) is much smaller than the width of the sec- 
ond (first) excited band, A2 (Ai): for 25 En (with 
uj z = 2tt x 400 Hz), V£/h » 100 Hz, Vf/h « 30 Hz, 
while A 2 /h w 16 kHz, and A1//1 « 1.6 kHz. Thus dipo- 
lar relaxation should be described by a Fermi golden rule. 
It is important though to take into account the coupling 
of the initial state to all energetically accessible states 
in the excited bands of the lattice, i.e. to all states 'J™ 
within an energy band qualitatively set by A,: each of 
the excited vibrational states (in the i th band) along the 
lattice axis can be considered as a continuum of states of 
width A^ and the coupling of the initial state to these 
many continua can be described as a sum of Fermi golden 
rules [la |. In the specific case where channel 2 is domi- 
nant this leads to: 

N = ^j2^ n o DN ( 4 ) 



Interestingly, the obtained ID dipolar relaxation rate 
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FIG. 4: (Color online) Heating in the lattice below threshold after 
40 ms, as a function (a) of the angle B between the magnetic field and 
the tube axis, (b) of the transverse ellipticity (uj x — uj y ) 2 / (lj x -\-cj y ) 2 of 
the lattice sites (where uJ x>y are the oscillation frequency of the lattice 
site along the two main axes x and y). Lines arc guides for the eye. 



not produced in our experiment. Using deeper lattices 
should strongly reduce tunneling and enable the coher- 
ent excitation (when V 2 ° > A 2 , at typically 75 Er) of 
such rotating states in the 2D lattice. Operating in 3D 
optical lattices in the Mott regime with two atoms per 
site would also provide a resonant character; as collisional 
de-excitation in 3D lattices is strongly reduced [f|, one 
may reach an intriguing situation where an insulating 
Mott state in the ground band and a superfluid state in 
the excited band are coupled by DDIs. 
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parameter fa depends on the trap ping parameter po, sim- 
ilarly to what we found in 2D [15] . In addition, as noticed 
above, fa depends on the width of the excited band: in 
that respect, resonant dipolar relaxation at threshold in 
the lattice is not one-dimensional. For 25 En, we find 
fa,th ~ 2 x 10~ 19 m 3 s _1 . To relate the measured heat- 
ing rate of 215 nK.ms -1 to a rate parameter, we average 
over the density of the collection of the ID gases (see 
[l§] for a similar calculation in 2D), and, assuming that 
dipolar relaxation occurs solely through channel 2, find 
fa, exp ~ (4.6 ± 1.4) x l(T 20 m 3 s" 1 . The slight disagree- 
ment between fajh and fa,exp can be due to the correla- 
tions in the ID quantum gases: in our situation the ID 
interaction parameter 7 [20j is close to one. 

Below threshold, we measure an extremely low heating 
rate, and we deduce /3 = (5±1.5) x 10~ 22 m 3 s _1 . This is 
typically three orders of magnitude smaller than for 3D 
degenerate quantum gases [15j |. As we observe almost 
no heating for up to 100 ms, this strong reduction en- 
ables us to reach a regime where the ID quantum Bose 
gas is metastable, despite its highly out-of-equilibrium 
magnetic state. This important result opens the way to 
stable spinor mixtures in ID. 

As shown in Fig. [4j both the cylindrical symmetry of 
the lattice sites and the proper alignment of the mag- 
netic field along this symmetry axis are needed to reach 
a regime where dipolar relaxation is strongly suppressed 
(at low B). The observed suppression below threshold is 
thus directly related to cylindrical symmetry and conser- 
vation of angular momentum. Indeed, dipolar relaxation 
above threshold should populate rotating states in the 
lattice, in the spirit of the EdH effect. However, due to 
tunneling in higher lattice bands the rotating states are 
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